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This invention was made with Government support under contract number 
DE-AC02-98CH10886, awarded by the U.S. Department of Energy. The 
Government has certain rights in the invention. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to positron emission tomography (PET) 
scanners, and more particularly to low-power, low-noise compact PET scanners for 
use in obtaining an input function from a portion of the human body, such as the 
10 wrist. 

Description of the Prior Art 

Positron emission tomography (PET) scanning is a diagnostic tool for non- 
invasive^ imaging living organisms. It remains essential to the investigation of 
chemical and functional processes in biochemistry, biology, physiology, anatomy, 
molecular biology, and pharmacology. While technique s, such as x-rays, computed 
tomography (CT), and magnetic resonance imaging (MRI) provide anatomical 
images, PET scanning provides insight into biochemical changes that generally occur 
long before a corresponding structural change is detectable by more traditional 
techniques. 

20 Positrons are positively charged electrons emitted by the nucleus of an 

unstable radioisotope. The radioisotope is unstable because it is positively charged 



and has too many protons. Release of the positron stabilizes the radioisotope by 
converting a proton into a neutron. For radioisotopes used in PET, the element 
formed from positron decay is stable. All radioisotopes used in PET decay by 
positron emission. The positron travels a small distance, which depends on its energy, 
before combining with an electron during a so-called "annihilation event". The 
annihilation event ultimately converts the mass of the positron and electron into two 
gamma rays that are emitted at 180° to each other along a so-called "line of 
coincidence". These gamma rays are readily detectable outside the human body. 

A small amount of the unstable radioisotope is typically administered to the 

10 patient by injection or inhalation, following which it circulates through the body. 
Scintillation crystals in a tomograph detect the gamma rays emitted by the 
radioisotope and convert them to light photons. The light photons are then converted 
to electrical impulses that are processed by the tomograph to determine the location of 
the annihilation event as being along the line of coincidence. 

Kinetic imaging or modeling measures the uptake of a tracer isotope over a 
period of time. The distribution of tracer isotopes may be used to represent regional 
blood flow and glucose metabolism. These studies often require catheterization to 
obtain discrete blood samples, which are analyzed for radioactivity and radioisotope 
metabolism. Unfortunately, the invasive withdrawal of blood is a significant 

20 discomfort to the patient, as well as a significant health risk for both the patient and 
hospital personnel through exposure to blood borne diseases and radioactive 
contamination. Therefore, direct arterial blood sampling is considered a health risk 
for both patients and health workers. 
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A wide range of quantitative PET studies using tracer kinetic modeling require 
accurately measured radiotracer concentrations in arterial blood as a function of time 
after injection, which is commonly referred to as an "arterial input function". To 
circumvent the health risks associated with direct arterial blood sampling, several 
approaches have been examined in an effort to non-invasively obtain an accurate 
arterial input function. While some approaches have focused primarily on the use of 
tomography, others have examined additional detector systems that generate a 
quantitative image-derived input function. 

Using the tomograph, studies have examined the possibility of obtaining an 

10 input function using large blood vessel imaging. However, this approach is limited in 
several respects. First, tomography exhibits a partial volume effect defined by spatial 
resolution. However, an artery large enough to provide reliable data may not be in the 
field of view. Second, time resolution may be determined by frame acquisition rates 
specified for a particular study. Although list mode acquisition capabilities reduce 
restrictions associated with slower acquisition rates, many scanners do not have this 
capability. Third, subject placement within the tomograph may affect the accuracy of 
the input function and obtaining reproducible positioning of the body is difficult. 

An alternative approach involves placing a radioactivity detector directly over 
a blood vessel or lung. The primary disadvantage of this approach is the substantial 

20 background associated with the surrounding tissue. This background must be 
subtracted to obtain the true input function. Since this approach is not based on 
coincidence counts, the signal may include a substantial amount of noise. 

Another alternative is to use a standardized input function, which is averaged 
across many subjects, or a modeled input function. In the latter method, the input 
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flinction is calculated from various physiological parameters. However, since the 
input function is very dependent on individual physiological states and procedural 
variables, such as differences in injection rates, these methods may lead to inaccurate 
results. Therefore, each of the techniques discussed above yield potential errors and 
there is a distinct need to determine accurate input functions by measuring the blood 
activity with little background from remaining portions of the body. 

Compact PET detectors require an efficient method of transmitting signals 
from the detectors to remote electronics for off-line processing. Conventional 
detectors communicate via independent data links, each of which is dedicated to a 
10 particular channel. However, since the majority of PET detectors include hundreds or 
even thousands of channels, this technique is too cumbersome for a compact PET 
detector. 

OBJECTS AND SUMMARY OF THE INVENTION 
It is an object of the present invention to provide a method and apparatus that 
reduce the discomfort of a patient, the risk associated with blood borne pathogens, 
and radioactive dispersion from the collection of blood samples typically required in 
kinetic modeling Positron Emission Tomography (PET) studies. 

It is another object of the present invention to provide a low-noise, low-power 
method and apparatus that are capable of non-invasively, selectively, and accurately 
20 measuring arterial radioactivity. 

It is yet another object of the present invention to provide a method and 
apparatus that maximize detector efficiency and spatial resolution by separating 
arterial information from venous and surrounding tissue information during PET 
scanning. 
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It is a further object of the present invention to provide a method and 
apparatus that utilize improved timing resolution and selective shielding techniques to 
effectively reduce random coincidence during PET scanning. 

A method of serially transferring annihilation information in a compact 
positron emission tomography (PET) scanner used to obtain an input function from at 
least a portion of a human body in accordance with one form of the present invention, 
which incorporates some of the preferred features, includes the steps of inputting a 
first time pulse on a first channel, and generating a first time signal representing the 
time-of-occurrence of the first time pulse and a first annihilation event. The method 
10 also includes generating a first address signal including a first address representing the 
first channel, generating a first channel signal including the first time signal and the 
first address signal, and outputting the first channel signal serially. 

The first time pulse is asynchronous to a clock signal and the first address 
signal is synchronous to the clock signal. The method may also include the steps of 
generating a composite signal including the first channel signal and a similarly 
generated second channel signal concerning a second annihilation event, and 
outputting the composite signal serially. 

An apparatus to serially transfer annihilation information in a compact 
positron emission tomography (PET) scanner used to obtain an input function from at 
20 least a portion of a human body in accordance with one form of the present invention, 
which incorporates some of the preferred features, includes a first time signal 
generator, a first address signal generator, and a first channel signal generator. The 
first time signal generator inputs a first time pulse on a first channel, which includes a 
position representing a time-of-occurrence of a first annihilation event. The first time 
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signal generator generates a first time signal representing the time-of-occurrence of 
the first time pulse and the first annihilation event. The first time pulse is 
asynchronous to a clock signal. 

The first address signal generator generates a first address, which represents 
the first channel, and a first address signal, which includes the first address. The first 
address signal is synchronous to the clock signal. The first channel signal generator 
generates a first channel signal including the first time signal and the first address 
signal, and outputs the first channel signal serially. 

A compact positron emission tomography (PET) scanner used to obtain an 
10 input function from at least a portion of a human body, which incorporates some of 
the preferred features, includes a scintillation array, detection array, front-end array, 
and serial encoder. The scintillation array includes a plurality of crystals and outputs 
photons in response to receiving gamma radiation from an annihilation event. The 
detection array includes a plurality of detectors and outputs a detection signal in 
response to detecting the photon. 

The front-end array includes a plurality of front ends and outputs a time pulse 
in response to receiving the detection signal. The serial encoder includes a time 
signal generator, address signal generator, and channel signal generator. 

These and other objects, features, and advantages of this invention will 
20 become apparent from the following detailed description of illustrative embodiments 
thereof, which is to be read in connection with the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 is a block diagram of a Positron Emission Tomography (PET) 
system, which includes a wrist detector formed in accordance with the present 
invention. 

Figure 2a is a pictorial view of a preferred embodiment of the wrist detector 
formed in accordance with the present invention. 

Figure 2b is a pictorial view of an alternative embodiment of the wrist detector 
formed in accordance with the present invention. 

Figure 3 is a pictorial view of an LSO crystal array and an APD detector array 
10 used for measurement of an input function in the wrist detector formed in accordance 
with the present invention. 

Figure 4 is a top outline view of the alternative embodiment of the wrist 
detector shown in Figure 2b. 

Figure 5 is a perspective outline view of an alternative embodiment of the 
wrist detector formed in accordance with the present invention. 

Figure 6 is a block diagram of a conventional technique for transferring 
information from PET detectors for remote processing. 

Figure 7 is a block diagram of a technique in accordance with the present 
invention for serially transferring information in accordance with the present 
20 invention. 

Figure 8 is a timing diagram showing signals associated with the serial transfer 
of information from PET detectors in accordance with the present invention. 
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Figure 9 is a graph of inefficiency or the probability that an event will be 
blocked by a busy link condition in the wrist detector formed in accordance with the 
present invention. 

Figure 10 is a block diagram of a preferred embodiment of the circuitry 
associated with detector channels in the wrist detector formed in accordance with the 
present invention. 

Figure 1 1 is a partial schematic diagram of a preferred embodiment of front- 
end circuitry and a block diagram of the remaining circuitry associated with detector 
channels in the wrist detector formed in accordance with the present invention. 
10 Figure 12 is a schematic diagram of an alternative embodiment of serial 

encoding circuitry associated with detector channels in the wrist detector formed in 
accordance with the present invention. 

Figure 12a is a schematic diagram of a preferred embodiment of the serial 
encoding circuitry shown in Figure 12 implemented in an Application Specific 
Integrated Circuit (ASIC). 

Figures 13-22 are schematic diagrams of an alternative embodiment of an 
ASIC implementation of the serial encoding circuitry formed in accordance with the 
present invention. 

Figure 23 is a block diagram of an experimental measurement setup for use in 
20 collecting sample data in accordance with the present invention. 

Figure 24 is a wrist phantom for use in collecting sample data in the setup 
shown in Figure 23. 

Figure 24a is a graph of an input function obtained by the wrist detector 
formed in accordance with the present invention. 
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Figure 25 is an arterial blood curve obtained using an automatic blood- 
sampling device. 

Figure 26 is a rectilinear image of the wrist phantom shown in Figure 24. 

Figure 26a is a block diagram of a test setup used to measure radioisotope 
activity between two detector blocks. 

Figure 27 is a graph representing the timing resolution or Root Mean Square 
(RMS) zero crossing jitter from the wrist detector formed in accordance with the 
present invention. 

Figure 28 is a distribution of cross-talk signals for a selected crystal in a 
10 scintillation array. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
A quantitative positron emission tomography (PET) imaging experiment 
generally requires inserting an arterial catheter, taking blood samples, and analyzing 
these samples for radioactivity concentration. The method and apparatus formed in 
accordance with the present invention measures the radioactivity present in the blood 
without the use of a catheter to physically remove the blood from the artery and 
mathematically corrects the results of this measurement to accurately represent an 
arterial input function. 

The wrist detector system formed in accordance with the present invention 
20 provides an accurate arterial input function in the presence of background activity 
from the body. Figure 1 is a block diagram showing the wrist detector 10 used in 
conjunction with a PET scanning system to determine the timing and location of 
annihilation events. 
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The system preferably includes a coincidence processor 12 that receives data 
from the wrist detector 10 through a Time-to-Digital Converter (TDC) 19 and a wrist 
detector cable 14. The TDC 19 preferably detects a timing edge with corresponding 
channel and block address information from the wrist detector 10 -corresponding to a 
particular annihilation event. This information is then stored by the coincidence 
processor 12. Power is preferably provided to the wrist detector 10 through the wrist 
detector cable 14 by a Low Voltage (LV) power supply 16 and a High Voltage (HV) 
power supply 18. The TDC 19 translates asynchronous timing pulses representing the 
time-of-occurrence of an annihilation event into digital information, such as a 

10 timestamp, which is used by the coincidence processor 12 to establish coincidence. 

Figure 2a is a pictorial view of a preferred embodiment of the wrist detector 
10, which includes two arrays of Lutetium Oxyorthosilicate (LSO) scintillation 
crystals 20 and two corresponding arrays of avalanche photodiode (APD) detectors 
22. An LSO crystal array 20 and an APD detector array 22 are also shown in Figure 
3. A conventional technique has been used with Bismuth Germanate (BGO) crystals 
and photomultipliers, which exhibits lower spatial resolution. However, such a device 
exhibits substantial inaccuracies as a result of large partial volume effects. LSO 
crystals have a detection efficiency similar to that of BGO, but have a faster 
scintillation decay interval, which is critical for reducing deadtime losses. 

20 In addition, the high stopping power of LSO crystals is ideal when using 

smaller crystal elements, which improves the spatial resolution of the detector. 
Spatial resolution is critical for imaging small arteries in the wrist. It is to be noted 
that although the description of the present invention is directed to a wrist detector, 
with modifications in size, placement, resolution, and other parameters, the detector 
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may alternatively be used to image additional portions of the human body, such as the 
head, neck, arm, and leg while remaining within the scope of the present invention. 

The spatial resolution of the LSO crystal array 20 approaches theoretical limits 
due to positron range, which is about 1 mm. The compact size of the APD detector 
array 22 is also an advantage when used in the limited space of a clinical scanner. As 
shown in Figure 3, the one-to-one correspondence between the elements of the LSO 
crystal array 20 and the APD detector array 22 is advantageous in identifying the 
point of gamma conversion. Thus, the compact size of each of the APD detectors 
combined with the potential for mounting multiple APD detectors in an array, provide 

10 a unique and compact multi-module detector that is not conventionally achievable 
with larger photomultiplier block detectors. 

As shown in Figure 2a, the wrist detector 10 formed in accordance with the 
present invention preferably includes two (2) LSO crystal arrays 20 and two (2) APD 
detector arrays 22, although any number of LSO arrays 20 and corresponding APD 
detector arrays, such as twelve (12) of each as shown in Figure 2b, is intended to be 
within the scope of the present invention. The wrist detector 10 shown in Figure 2b 
may be used to provide three-dimensional imaging of a portion of the human body 
including but not limited to the wrist. 

Each LSO crystal array 20 preferably includes thirty-two (32) LSO crystals 

20 configured as a 4 x 8 array. Each crystal is preferably about 2x2x10 mm 3 and is 
coupled to a corresponding APD detector element, an array of which is available from 
Hamamatsu Corporation, Bridgewater, New Jersey 08807 as part number S8550. 
The LSO crystal arrays 20 and APD detector arrays 22 are preferably mounted 
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together in a one-to-one configuration using a thin layer of silicone room temperature 
vulcanizing rubber adhesive or RTV. 

As shown in Figures 2a and 2b, the APD detector arrays 22 are preferably 
mounted on rigid Printed Circuit Boards (PCB) 24, which are coupled to each other 
by flex cable 26. Gamma radiation emitted by an annihilation event is first detected 
by crystals in the LSO crystal arrays 20, which translate the gamma radiation into 
photons. The photons are then detected by corresponding elements in the APD 
detector arrays 22, which translate the photons into an electrical signal. The electrical 
signal is coupled to electronic circuitry in a corresponding Application Specific 
10 Integrated Circuit (ASIC) 28, which is preferably mounted on a corresponding (PCB) 
24. Electronic signals from each of the ASICs 28 are preferably made externally 
available on a connector 30 through the wrist detector cable 14. 

Figure 4 is a top outline view of the wrist detector 10 shown in Figure 2b. A 
perspective outline view of an alternative embodiment of the wrist detector 10, which 
includes stacked pairs of LSO crystal arrays 20 and APD detector arrays 22 in each of 
twelve (12) detector assembly blocks 31, is shown in Figure 5. 

Time-of-occurrence information concerning an annihilation event and the 
corresponding address of the specific APD detector element that detected the 
annihilation event has conventionally been transferred to remote electronics, such as 
20 the TDC 19, for off-line processing, via separate lines. Each of these separate lines is 
dedicated to a single detector channel, as shown in Figure 6. The majority of PET 
scanning systems include hundreds or even thousands of these channels, which render 
the conventional method far too cumbersome for a compact imaging system. 
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The wrist detector formed in accordance with the present invention solves this 
problem by using a limited number of data links to transmit information concerning 
multiple annihilation event from a large number of channels (N) that share the link, as 
shown in Figure 7. The information required for each annihilation event is its time- 
of-occurrence, the address of the channel that detected the annihilation event, and the 
detector block corresponding to this channel. Information concerning the energy 
recorded in that channel may also be transmitted with the time and address 
information from the wrist detector in accordance with the present invention. 

The time-of-occurrence of the annihilation event is preferably represented by 
10 an asynchronous position of a leading edge of a data packet from each of the ASIC 
devices in the wrist detector 10. The address corresponding to the channel that 
detected the annihilation event is preferably serially encoded in the same data packet 
by a digital word having a length of log 2 (N) bits, where N is the number of channels 
in the detector block. Thus, for a detector block having thirty-two (32) channels, a 
minimum of five (5) bits would be required to encode the channel address. 
Information concerning the energy of the annihilation event is preferably encoded by 
the position of a second edge in the data packet. Therefore, edges representing time- 
of-occurrence and energy are preferably asynchronous with respect to a system clock 
while edges representing the channel address are preferably synchronous with respect 
20 to the system clock, as shown in Figure 8. 

Since the data link is shared between N channels, the total duration of the data 
packet is preferably much less than the average inter-arrival time of annihilation 
events in the N channels. This is preferably accomplished by increasing the 
frequency of the system clock. For instance, assuming no energy information is 
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transmitted and there is one guard bit, which has a duration of one period of the 
system clock, between the leading edge representing the time-of-occurrence and the 
beginning of the channel address, the maximum duration of the data packet is 
preferably represented by the following equation: 

^,=(log 2 (AO + 2)*7;, oa (1), 

where T clock is the period of the system clock. To prevent annihilation events from 
being blocked by a busy link condition the duration of the packet is preferably 
represented by the following equation: 

T pQCket {{\l{N*rate) (2), 

10 where rate refers to the average rate of annihilation events per channel. The busy link 
condition, when there is a conflict between the detection of two or more substantially 
simultaneous annihilation events, is preferably resolved by a priority encoder, which 
preferably neglects the event associated with the lower channel address. 

A graph showing the probability that an event will be blocked by a busy link 
condition, which is also referred to as inefficiency, is shown in Figure 9 for the case 
of N=32 channels and system clock frequencies of 50, 100, 200, and 500 MHz as 
indicated by lines 32, 34, 36, and 38, respectively. The abscissa of the graph is the 
singles rate for the entire scanner, which is assumed to have 384 total channels 
organized in twelve (12) blocks of thirty-two (32) channels. 

20 Figure 10 is a block diagram of a preferred embodiment of electronic circuitry 

corresponding to the channels of the wrist detector 10 formed in accordance with the 
present invention. The output of an element 40 of the APD detector array is 
preferably connected to a preamplifier 42, which inputs gain adjustment and 
calibration signals obtained from a digital-to-analog converter (DAC) 43. The output 
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of the preamplifier 42 is preferably shaped to be about a 70ns pulse by one or more 
filtering networks in a shaper network 44. The shaped pulse is preferably inputted to 
a Zero Crossing Detector (ZCD) or Constant Fraction Discriminator (CFD) circuit 46, 
which generates precise logic pulses in response to the shaped pulse input crossing a 
given threshold, as well as the peak of the shaped pulse. 

Detection of the peak of the shaped pulse preferably yields a timing pulse 
from the ZCD/CFD 46 that represents a time-of-occurrence of the corresponding 
annihilation event. The time between the first threshold crossing of the shaped pulse 
and the second threshold crossing of the shaped pulse indicates the energy of the 

10 annihilation event using a so-called "time-over-threshold method". The energy of the 
annihilation event is preferably represented by the position of an energy pulse from 
the ZCD/CFD 46. Since the energy measurement is being performed on the shaped 
pulse following the preamplifier 42 and shaper circuit 44, the value of the energy 
measurement is preferably calibrated to yield a more accurate result. 

A DAC 48 preferably provides a threshold signal and a threshold adjustment 
signal to the ZCD/CFD 46. The ZCD/CFD 46 is preferably based on two 
comparators. One comparator is used for arming the ZCD/CFD 46 by triggering on 
signals having the proper energy and independently setting the threshold signal for the 
detector channels. The threshold is preferably set for the detector channels through a 

20 serially loaded shift register (not shown). The shift registers for each of the detector 
block ASICs are preferably daisy-chained and share the same data, control, and clock 
signals. The shift registers preferably set the detector channel DACs 43, 48 to 
appropriate values. The remaining comparator in the ZCD/CFD 46 is preferably used 
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to determine the baseline crossing of the bipolar signal, which represents the energy 
of the annihilation event as described above. 

CFDs generally use a constant fraction or percentage of the input pulse to 
determine the timing of the output pulse relative to the input signal. This technique is 
not subject to jitter, which is typically caused by varying amplitudes or rise times of 
the inputs, such as in leading edge discriminators. The pulse output from the 
ZCD/CFD 46 preferably has a standardized amplitude and a preset duration. 

The output of the ZCD/CFD 46 is preferably applied to an edge/address 
encoder 50 that outputs a serialized edge/address signal 74, which includes the time- 
10 of-occurrence and the channel address corresponding to the detected annihilation 
event. The system clock signal 54 is preferably inputted to the edge/address/priority 
encoder 50 for synchronous timing. An analog monitor multiplexer 56 is preferably 
provided to enable external access to analog signals within the channel circuitry, such 
as the analog signals before and after the ZCD/CFD 46. 

Figure 1 1 is a schematic diagram of front-end circuitry in relation to 
subsequent functional blocks for each of the detector channels. Preamplifier and 
shaping network parameters are preferably optimized with respect to technological 
parameters and operating point characteristics of the wrist detector, such as leakage 
current and gain, to minimize the Equivalent Noise Charge (ENC). Simulations have 
20 predicted an ENC of about 700 electrons RMS at about a 70ns peaking time. A gain 
of 3.3 mV/fC (where Cf = 300fC) was set for the preamplifier with about 1 .3mW of 
power dissipation. 

The output of the APD element 40 is provided to the preamplifier 42 that 
preferably includes a parallel connection of a capacitor CI and a resistance Rl, which 
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are connected in parallel across an inverting terminal and an output terminal of an 
operational amplifier 58. The resistance Rl is preferably realized by a Field Effect 
Transistor (FET) and the preamplifier 42 is connected in series with the output of the 
APD element 40. 

The output of the preamplifier 42 is provided to the shaper circuit 44, which 
preferably includes a pole-zero cancellation network, a l st -order low pass filter, and a 
2 nd -order bandpass filter. An analysis was performed to determine the peaking time 
that optimizes timing resolution. Assuming that the LSO decay time is 40ns and 
considering the series noise, the time from the peak to the zero crossing of a bipolar 

10 Gaussian between 25ns to 90ns would lead to a theoretical optimum of about 700ps 
rms zero-crossing jitter, as indicated in Figure 27. Therefore, the peaking time is 
preferably set to about 70ns. Each transistor in the front-end circuitry is preferably 
optimized to minimize its electronic noise contribution. 

The pole-zero network preferably includes a parallel combination of one or 
more capacitors C2 connected in series and one or more resistances R2 connected in 
series. The number of capacitors C2 is preferably equal to the number of resistances 
R2, as well as being greater than two (2). The resistance R2 is preferably realized by 
an FET and the pole-zero network is connected in series with the output of the 
preamplifier 42. The pole-zero network is preferably used to compensate the reset 

20 transistor non-linearity, as well as to reduce the noise contribution from subsequent 
stages. 

The l st -order shaper network preferably includes a parallel combination of a 
capacitor C3 and a resistor R3 connected in parallel across an inverting terminal and 
an output terminal of an operational amplifier 60. A resistor R4 and a capacitor C4 



-18- 

are preferably connected in series between an output terminal of the operational 
amplifier 60 and an inverting terminal of an operational amplifier 62. A capacitor C5 
and a resistor R5 are preferably connected in parallel across the capacitor C4 and the 
resistor R5 is connected in parallel across the inverting terminal and the output of the 
operational amplifier 62, which is essentially the 2 nd -order shaper network. 

The non-inverting terminal of each of the operational amplifiers 58, 60, 62 is 
preferably connected to ground and an output terminal of operational amplifier 62 is 
provided to the ZCD/CFD 46. An N-channel edge/address/priority encoder circuit 64 
is shown in Figure 1 1, which preferably includes the edge/address/priority encoder 

10 circuits 50 shown in Figure 10 for N channels in a detector block. The serial encoder 
circuit 64 outputs a single serial output signal 52 to the TDC. Further details 
concerning the front-end electronics for the detector channels are provided in U.S. 
Patent No. 5,793,254, which is incorporated herein by reference. 

Figure 12 is a schematic diagram showing an alternative embodiment of the 
edge/address encoder 50 shown in Figure 10, which may be repeated for each of the 
detector channels. A CFD signal 66, which includes a pulse representing the time-of- 
occurrence of an annihilation event is applied to two (2) inverters HS5, HS6 for initial 
pulse shaping. The output of inverter HS6 is used to asynchronously clock a high 
level (Vddd) through a flip-flop HS1, which provides an edge signal 68 outputted 

20 from flip-flop HS 1 . The edge signal 68 is synchronously clocked through two 
additional flip-flops HS2, HS4 and the inverted output of flip-flop HS4 is used to 
reset flip-flop HS 1 . The inverted output of flip flop HS2 and the inverted output of 
flip flop HS1 are applied to the inputs of a nor gate HS30, the output of which is used 
as a load signal 70 for a shift register (HS34, HS33, HS36, HS37, HS38). 
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In the embodiment shown in Figure 12 the channel address is provided as a 5- 
bit digital quantity. Each of the A inputs of multiplexed flip flops HS34, HS33, 
HS36, HS37, HS38 is connected to a high level (Vddd) or a low level (Vssd) to 
represent corresponding bits of the channel address. For instance, if the channel 
address is defined by bits A0-A4 (AO being the least significant bit) then a binary 
channel address of 1 1 101 is represented by connecting the A input of flip flop HS37 
to a low level and connecting the A inputs of the remaining flip flops HS34, HS33, 
HS36, HS38 to a high level, as shown in Figure 12. 

The load signal 70 is an active high pulse having a duration of about one 

10 period of the system clock signal 54. When the load signal 70 is high, the A inputs of 
each of the multiplexed flip flops HS34, HS33, HS36, HS37, HS38 are loaded into the 
corresponding flip flop. When the load signal 70 is low, the flip-flops HS34, HS33, 
HS36, HS37, HS38 function as a 5-bit serial shift register, the output of which is a 
channel address signal 72. The channel address signal 72 is combined with the edge 
signal 68 using a nor gate HS35 and inverters HS39, HS40, HS41 to yield an 
edge/address signal 74. 

A timing diagram of these signals is provided in Figure 8, which shows that 
the edge/address signal 74 is a combination of asynchronous analog information 
representing the time-of-occurrence of an annihilation event by the position of an 

20 edge, as well as synchronous digital information representing the address of the 
particular channel that detected the annihilation event. The combination of this 
information permits a single serial data link to be shared by each of the channels in a 
detector block or each of the channels in the wrist detector formed in accordance with 
the present invention. 
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Figures 13-22 are schematic diagrams of an alternative embodiment of the N- 
channel serial encoder 64 shown in Figure 1 1, in which N=32. These diagrams were 
prepared in contemplation of an ASIC implementation of the encoder 64, which 
essentially incorporates the encoder circuitry shown in Figure 12 for each of the 32 
channels, as well as circuitry to combine the outputs of each of the channel encoder 
circuits into a single serial output signal 52 shown in Figure 1 1 and a priority encoder 
to resolve busy link conditions. The ASIC is preferably designed to have minimal 
power dissipation so as not to affect the gain of the APD detector array, which is 
sensitive to temperature. 

10 Figure 12a is a schematic diagram of a preferred embodiment, which is 

implemented in an ASIC, of the serial encoding circuitry shown in Figure 12. The 
preferred embodiment is similar to that shown in Figure 12, except that a single shift 
register 67 is used to serialize the addresses for each of the channels in the ASIC. The 
channel address A0-A4 is preferably generated by a 32-to-5 priority encoder, such as 
an encoder 39 shown in Figure 13. 

An asynchronous trigger signal 71, which represents the time-of-occurrence of 
an annihilation event, is preferably applied to a time signal generator 73. A 
serializing signal 75, which is a synchronous derivative of the trigger signal 71, 
propagates through the shift register 67 to generate enabling signals Select0-Select4. 

20 The enabling signals Selecte0-Select4 gate successive bits of the channel address 
through combinatorial logic 77 to a flip flop 79 that outputs the serialized channel 
address signal 72. 

In an experimental measurement setup shown in Figure 23, which is intended 
to represent an alternative embodiment of the present invention, signals from the APD 
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detector array 22 are preferably inputted into a charge sensitive preamplifier 76 and 
then into shaping amplifiers 78 with about a 70 ns peaking time. The output of the 
shaping amplifiers 78 are then preferably divided between constant fraction 
discriminators (CFD) 80 with about a 50mV threshold and Computer Automated 
Measurement and Control Fast Encoding and Readout analog-to-digital converters 
(CAMAC FERA ADC) 82 that digitize these signals. A Nuclear Instrumentation 
Methods (NIM) logic unit or coincidence logic/gate generator 84 preferably selects 
coincidence events and generates a gate for the ADC's 82. Data are preferably 
transferred from the ADCs 82 and analyzed in an external system 86. 

10 Sensitivity of the APD detector array 22 is preferably determined using a 

known level of radioactivity in a 2 mm line source. An estimate of this sensitivity is 
also preferably calculated for the volume of radioactivity in the field of view of the 
detectors 22, which is about 0.05cm 3 , and compared to the measured value. The 
measured value is about O.lOHz/nCi. The count rate is typically low with this 
sensitivity, but is preferably increased by integrating over several seconds. It should 
be noted that the described sensitivity relates to a single pair of detectors that are 
about 10 mm deep. The sensitivity is preferably enhanced by increasing the quantity 
of detectors in the array and/or the depth of the LSO crystal arrays 20. 

An energy resolution study is preferably performed using the APD detector 

20 array 22. The energy resolution is preferably sufficient to clearly delineate a 5 1 1 keV 
signal from the background. LSO crystal arrays 20, which are available from Proteus, 
Inc., Chagrin Falls, Ohio 44022 and CTI, Inc., Knoxville, Tennessee 37932, are 
preferably used in the experiment. 
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LSO crystal arrays 20 available from CTI provide a signal of about 2500 
primary photoelectrons per MeV and an energy resolution of about 23% Full Width 
Half Maximum (FWHM) for 5 1 1 keV gamma rays. Two configurations were tested 
using the LSO crystal arrays 20 available from Proteus, each having a reflective 
barrier between the elements, which has excellent reflective properties and improves 
energy resolution. In one configuration, the reflector was bonded to the crystal 
surfaces, and in the other, it was not. 

The bonded array provides about 2300 primary photoelectrons per MeV and 
an improved energy resolution of about 17%, while the unbonded surface provides 
10 about 2600 primary photoelectrons per MeV and an energy resolution of about 13%. 
The APD array 22 provides an average gain of about 50, which results in a signal-to- 
noise ratio of at least 50:1. Differences in the gain of individual channels caused a 
shift in the associated pulse height spectra. These differences are preferably 
compensated for by shifting the peak position of the photopeak in each channel to the 
average peak position of the thirty-two (32) channels in each array. This makes it 
more convenient to set a threshold such that only photopeak events are used for image 
reconstruction. 

In order to determine the effect of attenuation on the sensitivity of the detector 
array 22, a Lucite® block 88, as shown in Figure 24, is preferably used to simulate a 
20 human wrist. The block 26 preferably includes one or more holes in the center large 
enough to accommodate Tygon® tubing to simulate arteries. The block 88 is 
preferably about 4.5 cm in thickness, which is consistent with the average thickness of 
the human wrist. The tubing 90 preferably runs through the block 88, which is placed 
between the detector arrays 12, as shown in Figure 23. 
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A series of simulated input functions is preferably acquired. An example of an 
input function is shown in the graph of Figure 24a. The block 88 is then removed and 
similar data is preferably reacquired. Experimental results show that the only 
difference between the two configurations was the block 88 and attenuation effects 
from the block 88 produced a 27% reduction in the sensitivity of the measurements. 

The wrist phantom block 88 is preferably fitted with tubing 90 chosen to 
closely parallel the diameter of both arteries. A solution of F-18 fluoride is preferably 
pumped through the tubing 90 while data is acquired using the system described 
above. Using this phantom, a sensitivity of about 0.04 cps/nCi/cc was obtained. 
10 Simulations are preferably performed using estimates of the noise established 

from typical count rates associated with the dynamic phantom experiments described 
above. The purpose of these simulations was to determine the magnitude of the error 
introduced by using data with a similar uncertainty to that which would be expected 
from a single pair of crystal arrays used in the wrist detector formed in accordance 
with the present invention. 

In these experiments, an actual arterial blood curve, which was obtained using 
an automatic blood-sampling device and shown as a solid line 92 in Figure 25, was 
used to calculate global glucose metabolic values using the Sokoloff model. Noise 
based on the expected uncertainties determined in the experimental measurement of 
20 the count rates is then preferably added to the blood curve using a random noise 
generator, which is represented as error bars 94 in Figure 25. The simulated blood 
curves are preferably used to calculate metabolic values in the presence of noise 
introduced by the wrist detector. 
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Ten input function curves were generated based on the expected signal to 
noise ratio for the wrist detector. The values obtained for whole brain metabolism 
using these data showed about a 3% uncertainty in the glucose metabolic rate. This 
fits well within the 8% uncertainty typical of clinical PET studies. The error bars 94 
may be reduced by a factor of two by increasing the detectors by a factor of four. 

A planar image of the line source is preferably used to estimate the spatial 
resolution of the block detector system. A rectilinear image of the wrist phantom is 
shown in Figure 26. Radioactivity is preferably decayed in the tubing to simulate a 
long acquisition image. The image is preferably used to identify regions of interest in 
10 the wrist that can then be used to generate the input function. 

One of the most difficult problems in measuring arterial radioactivity is its 
separation from venous circulation. A typical adult radial artery is about 2 to 3 mm in 
diameter. The 2 mm resolution of the directly coupled detectors is preferably used to 
separate an arterial signal from a venous signal. The range within which the artery is 
distinguishable from the vein is preferably defined using the wrist phantom. Such an 
experiment is preferably performed by pumping water containing a positron emitting 
radionuclide, such as fluorine-18 or carbon-1 1, through silicone tubing that has been 
placed in the channels of the wrist phantom, as shown in Figure 26a. The 
radioactivity is then preferably diluted in a reservoir and returned through channels 
20 representing the veins. In this way, the arteries may be distinguished from the veins. 

In a clinical setting, detector acquisition generally begins with radiotracer 
injection. As the radioactivity moves through the arteries in the wrist, a region of 
interest is preferably defined using a two-dimensional image obtained from the 
arteries alone. Once this image has been obtained, the pixels corresponding to the 



-25- 

image are preferably included in the signal for the arterial input function. The 
experiments involving the wrist phantom provide information concerning the spatial 
resolution required to obtain data exclusively from the artery while minimizing 
interference from the vein. It is also anticipated that the venous signal may be 
subtracted from the mixed signal, thereby improving the accuracy of the arterial input 
function. 

Several approaches may be used to maximize the detected signal. The first is 
to determine the count rate from the wrist phantom using radioactivity levels in the 
tubing that are similar to levels of radioactivity present during a typical PET scan. 

10 The wrist phantom is preferably used to measure count rate and determine the 

optimum geometry that can be used to maximize signal-to-noise ratio. This yields an 
optimal position of the detectors so that the signal from the artery is high while the 
signal from the vein is low. The optimal position of the two detectors is preferably 
determined for both the radial and ulnar arteries. 

In an alternative embodiment, the wrist phantom preferably includes one pair 
of detectors on the ulnar artery and an additional pair of detectors on the radial artery. 
Yet another pair of detectors running along the length of both the radial and the ulnar 
arteries may also be used. This provides a sum of four detector pairs, each of which 
preferably counts in coincidence. The arterial input function is then preferably 

20 determined using the sum of the artery signal from the radial artery and the ulner 
artery. This approach maximizes the signal while preserving resolution to separate 
the arterial signal from the venous signal. 

Provided the geometry can be fixed and the attenuation generated by the wrist 
can be determined, quantitative information necessary to obtain a useful input 
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function is preferably available directly from the detector. The arterial input function 
may alternatively be normalized to venous samples obtained at a time when arterial 
and venous blood are in equilibrium with respect to radioactivity concentrations. 

Noise in the electronics and readout system of the wrist detector represents a 
major limitation to timing resolution. Therefore, the Application Specific Integrated 
Circuit (ASIC) shown in Figures 2a, 2b, 4, and 13-22 is preferably substituted for 
discrete devices in the wrist detector. The ASIC preferably uses 0. 1 8 \un 
Complementary Metal Oxide Semiconductor (CMOS) technology, as described in P. 
van Zant, Microchip Fabrication, 3 rd Edition, pp. 99-1 18 (1997), which is 
10 incorporated herein by reference, and is about 1.5x4.2mm 2 . 

Substitution of the ASIC significantly improves timing resolution. Current 
experimental data indicate a resolution of about 5 nanoseconds FWHM for the 
detector without the ASIC. Figure 27 is a graph representing the timing resolution or 
Root Mean Square (RMS) zero crossing jitter from the wrist detector with the ASIC 
formed in accordance with the present invention. 

As indicated in Figure 27, for about 1000-2000 photoelectrons produced with 

a gain of 50 in the APD detector array and a peaking time of about 70 nanoseconds, 

the timing resolution is about 1 nanosecond. Therefore, with a coincidence timing 
♦ 

resolution of 1 to 2 nanoseconds, much of the interference generated from random 
20 noise during acquisition of the input function can be eliminated. 

The light output of the detector array is preferably measured by coupling each 
array to the APD detector array. A 22 Na source is preferably positioned in front of the 
detector array and an energy pulse height spectrum is recorded for each channel. The 
peak position is then preferably measured, gain-corrected, and averaged for each LSO 
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crystal array. The unbonded reflector shows a higher light output than a bonded array 
having the same length and manufacturer. 

A reflective mask is preferably positioned on the end of the crystal array that 
is in contact with the APD. The light output of the detector is increased by allowing 
light only through an aperture that exactly matches the sensitive area of the APD 
detector array. A mask cut from reflective coating (available from 3M Corporation, 
St. Paul, Minnesota 55144) is preferably used to enhance light gain in the crystal 
arrays, which improves both energy and timing resolution. 

As shown in Figure 10, the ASIC preferably includes access to analog signals, 
10 which are preferably brought off the chip and used externally to measure energy. As 
shown in Figures 10 and 1 1, the preferred embodiment of the wrist detector includes a 
zero crossing discriminator (ZCD) or constant fraction discriminator (CFD) per 
channel and the data is transferred over a serial link from the ASIC on the wrist 
detector to a TDC near the patient. The TDC preferably adds an absolute time stamp 
and sends the time-of-occurrence, channel address, and block address information to a 
remote coincidence processor 12, as shown in Figure 1. A Versa Module Europa 
(VME) bus system is preferably used to increase the data transfer capability from the 
detectors to the remote coincidence processor 12. The VME bus preferably includes a 
TTL-based backplane which, although the system is asynchronous, sets the data 
20 transfer speed to about 20 Mbytes per second. 

A single incident gamma ray that causes signals in more than one readout 
channel may generate cross talk between crystals. Cross talk may reduce spatial 
energy, and time resolution, both of which are critical. Cross talk is typically caused 
by Compton scattering from one crystal to another or the escape of photoelectrons 
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from the primary crystal following absorption. In addition, for practical reasons, 
optical isolation of the scintillation photons within a single crystal is never perfect. 
Each of these effects may increase as crystal size is reduced. 

The degree of cross talk has particularly important implications in the design 
of front-end electronics. If all of the gamma-ray signal is contained within one 
crystal, there is theoretically no need for energy digitization since an energy window 
is enforced via simple discriminator thresholds for each crystal. 

One method that is preferably used to overcome cross talk is to sum the energy 
in the nearest neighbor crystals to recover the complete energy deposited by the 
10 photon. This method requires analog information to be collected concerning each 
detector. Simulations have been performed in accordance with this technique and the 
results are shown in Figure 28, which depicts the distribution of cross-talk signals for 
a selected crystal. A quantitative measurement of cross talk is the increase in signal 
produced by the inclusion of adjacent crystals. 

It is to be noted that there is less cross-talk between crystals separated by a 
continuous reflector sheet (along the long axis) than between crystals separated by 
narrow, cut reflector strips (along the short axis). There is about a 25% increase in 
coincidence sensitivity when using the combined signals, which is preferably 
achieved by adding the analog pulse height information to the readout chip. 
20 Improvements in the placement of reflective sheets also decrease the amount of cross 
talk. 

Concerning detector sensitivity, simulations performed using the expected 
signal-to-noise ratio generated from a single set of crystals provides about a 3% 
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average error in the glucose metabolic rate. Increased sensitivity and reduced 
uncertainty are preferably achieved by using additional and/or deeper crystals. 

Concerning background radioactivity from the nearby body and surrounding 
wrist tissue, specifically designed shielding combined with the high resolution of the 
wrist detector formed in accordance with the present invention selectively reduces or 
compensate for this interference. In addition, fast coincidence timing resolution 
minimizes these background effects. 

Concerning interference from venous blood, which is in close anatomic 
proximity to the arteries from which the input function is obtained, the high resolution 
1 0 of the detectors combined with simultaneous measurement of venous radioactivity 
and subtraction of this parameter from the arterial signal enhance the ability to 
accurately measure the arterial input function. 

Thus, the method and apparatus formed in accordance with the present 
invention may be used to non-invasively, selectively, and accurately measure arterial 
radioactivity curves for PET radiotracers with low-power requirements and low noise. 
The method and apparatus is able to generate an arterial input function that greatly 
reduces biohazards associated with blood withdrawal, as well as discomfort 
associated with arterial cannulation, which is useful in quantitative PET tracer studies 
in clinical research and diagnosis. The method and apparatus also maximize detector 
20 efficiency and spatial resolution by separating arterial information from venous and 
surrounding tissue information. 

Although illustrative embodiments of the present invention have been 
described herein with reference to the accompanying drawing, it is to be understood 
that the invention is not limited to those precise embodiments, and that various other 
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changes and modifications may be effected therein by one skilled in the art without 
departing from the scope or spirit of the invention. 



